Analytical model for electromagnetic cascades in rotating electric field 
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Electromagnetic cascades attract a lot of attention as an important QED effect that will reveal 
itself in various electromagnetic field configurations at ultrahigh intensities. We study cascade dy- 
namics in rotating electric field analytically and numerically. The kinetic equations for the electron- 
positron plasma and gamma-quanta are formulated. The scaling laws are derived and analyzed. For 
the cascades arising far above the threshold the dependence of the cascade parameters on the field 
frequency is derived. The spectra of high-energy cascade particles are calculated. The analytical 
results are verified by numerical simulations. 
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I. INTRODUCTION 



Quantum electrodynamics (QED) effects in a strong 
electromagnetic field attract a lot of attention P, [1] • In 
the last decade a new surge of interest in QED effects 
in superstrong electromagnetic field arises thanks to sub- 
stantial progress in laser technologies Theoreti- 
cians are considering such effects as vacuum polarization 
in laser field [6|], electron-positron (e+e~) pair produc- 
tion from vacuum in the field of colliding laser pulses 
and in combined Coulomb and strong laser field [sj etc. 

Up to now the electromagnetic fields are not so in- 
tense in laboratory to observe strong-field QED effects 
directly. However, the interaction of ultra-high energy 
electron beams with electromagnetic fields can be used to 
study the effects as in the frame of the relativistic electron 
the electromagnetic field can be very strong. The spec- 
tra of the scattered electrons and photons arising from 
the nonlinear Compton scattering of intense laser pulse 
by 46.6 GeV electron beam have been measured exper- 
imentally [13]. The e'^e pair production has been 
also observed. The strong-field QED effects like channel- 
ing radiation and e~^e~ pair production at the interac- 
tion of relativistic electron beam with crystalline fields 
has been also studied experimentally and theoretically 
[HI, [13 • ^ number of theoretical and numerical mod- 
els for QED effects accompanying an interaction of rel- 
ativistic electrons with strong electromagnetic fields has 
been developed. The evolution of the distribution func- 
tions of electrons, positrons and hard photons in QED- 
strong pulsed laser fields interacting with electron beams 
were found numerically [ISj] . In Ref . |14i | multiphoton tri- 
dent pair production in strong laser field was studied in a 
manifestly nonperturbative domain in the periodic plane 
wave approximation. In addition, the energy spectrum of 
positrons produced in collision of strong laser pulse with 
high-energy electron beam was calculated. The trident 
pair production amplitude in a strong laser background 



was also calculated in Ref. [T5|. As opposed to fl^, the 
finite pulse durations were allowed, while the laser fields 
were still treating nonperturbatively in strong-field QED. 
The individual contributions of the one-step and two-step 
processes were explicitely identified. 

Another effect that can be very important at high in- 
tensities is the generation of electromagnetic cascades 
[ll,[l3- The cascades develop as follows. A slow electron 
can be accelerated in superstrong laser fields up to very 
high energy. Then it can produce a high-energy photon 
7 by Compton scattering with n laser photons involved 
in the process: 



-|- nuj — ^ e' -f 7. 



(1) 



The resulting photon creates a pair in a photon- 
multiphoton collision: 



7- 



(2) 
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This reaction can be considered as the strong-field gener- 
alization of the Breit- Wheeler process [3, [3 • The pro- 
duced electron and positron as well as the initial (seed) 
electron are accelerated by the field and become capa- 
ble to produce new high-energy photons by the process 
dl]). Hence, the number of the particles grows rapidly 
due to the deposition of the laser energy into the par- 
ticle acceleration. As a result, the relativistic electron- 
positron-gamma plasma is produced. Besides pair cre- 
ation by consecutive processes (jlj and ([2]), the particle 
multiplication can occur by the process of Bethe-Heitler 
type (or trident process) [I, [11]: e + nuj ~^ e'e+e^. How- 
ever, the threshold intensity for cascade development is 
about / ^ 10^^ W/cm^, and for modern laser system the 
characteristic wavelength is about Xi 1 /im. At such 
parameters the overall probability of processes (HI) and 
(21) are much higher than the probability of the Bethe- 
Heitler process. 

The seed e+e^ pair can be created as a result of the in- 
stabihty of QED vacuum in strong electromagnetic field 
[3, [mil- The probability of the vacuum pair creation 
increases in the field of two colliding laser pulses 123, 23] 
or in the focused laser pulse [l^. However, the threshold 
intensity for the pair creation by colliding and focused 
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laser pulses with A; ~ 1 /im is very high and can be up 
to 2.3 X 10^^ W/cm^ for the pulse duration 10 fs and cir- 
cular polarization 26]. The use of high-energy photons 
as a seed for the cascade is more preferable than that of 
free electrons and positrons because the ponderomotive 
force of the laser pulses pushes the charged particles out 
from the high-intensity region. The seed particles can 
be also produced by atom ionization j27,] or pair produc- 
tion in the superposition of laser field and the Coulomb 
field of heavy nuclei Q. When the number of the cre- 
ated cascade particles becomes large the plasma effects 
are important and the significant amount of the laser 
energy can be absorbed due to production and heating 
of the electron-positron pair plasma. This can limit the 
attainable intensity of high-power lasers jTJ '28'] . 

The dimensionless parameters characterizing QED ef- 
fects can be introduced [1^, [13] : 



Xe 



eh 



(E + 



mc 
Vc,B])_ 



Ec 



-7, 



(3) 
(4) 



where is the 4 



vector of the field potential, E and B 
are the electric and magnetic field strength, respectively, 
F^v is the field-strength tensor, is the particle four- 
momentum, 7 is the particle Lorentz factor (7 = hio /mc^ 
for the photon), v is the particle velocity, E^r = rri^C' /eh 
is the QED critical field strength, e > is the value of 
electron charge, m is the electron mass, c is the speed 
of light, h is the Planck constant, uj is the photon fre- 
quency, LOi is the laser frequency, symbol _L denotes the 
component which is perpendicular to the particle veloc- 
ity. The first parameter oq characterizes the laser filed 
strength. If the laser field is week ao <C 1 then mostly 
n = 1 in Eqs. ([IJ and ©. In the opposite limit oq ^ 1 
the electron absorbs a large number of photons and the 
radiation spectrum becomes synchrotron- like (sOj- The 
second parameter Xe characterizes the regime of the pho- 
ton emission. In the limit Xe ^ 1 the radiation can be 
described classically (except the case when the number 
of the emitted photons is small, see Ref. [3l|) while in the 
opposite limit the recoil imposed on the electron by the 
emitted photon is substantial and the quantum approach 
should be used. Similarly, pair production is significant 
in the limit X7 ^ 1 while it is exponentially suppressed 
in the classical limit X7 ^ 1- 

At high intensity the so-called radiation formation 
length [201 becomes much smaller than the laser wave- 
length and the characteristic path passed by the elec- 
tron between two consecutive photon emissions. This 
regime is often referred as quasistatic or tunneling and 
can be treated in the framework of quasiclassical ap- 
proach [Hi, i29i]. The Compton scattering is treated just 
as the transition of the electron from one level to another 
in external constant and homogeneous electromagnetic 
field that is accompanied by photon emission. The pair 
production is considered similarly. In addition, the quasi- 
classical approach implies that the particles moves along 



the classical trajectories between the instances of pho- 
ton emission or pair production that strongly simplifies 
calculations. 

Quantum dynamics of the particles and their interac- 
tion are complex phenomena and a number of numerical 
models is proposed. There are numerical models based 
on the first principles like solving Dirac's equation and 
calculating the cross-sections of QED processes [32i]. To 
simulate collective dynamics of the particles the numer- 
ical solution of the kinetic equations can be used 
Another probably more numerically efficient approach is 
provided by Monte Carlo methods. A number of sophisti- 
cated models like well-known Geant4 toolkit for the sim- 
ulation of the passage of particles through matter [sl] 
are based on Monte Carlo methods. Monte Carlo calcu- 
lations are also used to analyze the cascade development 
in the strong electromagnetic fields (s^, [11]. Recently 
particle-in-cell methods have been also used to simulate 
the effect of the self-generated electromagnetic fields of 
the produced e^e~ pairs on the cascade development 

In this paper we examine and solve analytically the 
kinetic equations for the particle distribution functions 
in the electromagnetic cascade in the homogeneous ro- 
tating electric field. This is one of the most simple field 
configuration where electromagnetic cascades arise. Be- 
sides, the rotating electric field can be found in the B- 
node of the circularly polarized standing wave formed 
by two co-propagating circularly polarized laser pulses. 
The plane _B = is the most favorable for the cascade 
development. However, for the simplicity we consider 
homogeneous electric field. In order to analyze the par- 
ticle distribution we use the following assumptions: (i) 
we use the quasiclassical synchrotron formulae for the 
probabilities of photon emission and pair production, and 
neglect the angles < I/7; (ii) the probability rates are 
averaged over the polarization of the particles; (iii) the 
self-generated electromagnetic fields of the particles pro- 
duced in the cascade are neglected. The last assumption 
is valid at the initial stages of the cascade development 
when the density of the self-generated electron-positron 
pair plasma is not high [2^. Under the assumption that 
the number of the cascade particles doubles many times 
on the period of the field rotation the analytical solutions 
for the distribution functions are obtained for the region 
in the phase space 7 ^ (7) and x ^ 1: where (7) is the 
average relativistic factor of the cascade particles. The 
obtained solution is used to calculate the energy spectra 
of the electrons, the positrons and the photons. The an- 
alytical results are compared with the results of the nu- 
merical simulations performed using the particle-in-cell 
Monte Carlo technique described in Ref. [28] . 

The paper is organized as follows. In Sec. II we con- 
sider the electron, positron and photon motion in the ro- 
tating electric field. In Sec. Ill the cascade kinetic equa- 
tions for the particle distribution functions are discussed. 
The cascade scalings are derived and discussed in Sec. 
IV and Sec. V. The distribution functions of the cascade 
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particles are derived in Sec. VI. In Sec. VII the obtained 
results are verified by the numerical simulations and the 
asymptotic formulae for the particle spectra are given in 
Sec. VIII. Finally, in Sec. IX summary and discussion are 
given. 



II. PARTICLE MOTION IN THE ABSENCE OF 
PHOTON EMISSION AND PAIR PRODUCTION 

First we study particle motion in the rotating electric 
field in the framework of the classical approach, neglect- 
ing photon emission and e"'"e~ pair production. The field 
components are the following: 



= Eocos{t + 0o), 
Ey = Eosm{t + ^o), 

that leads to the following equations of motion 



dpx 
dt 
dpy 
lit 



-COs(t + (/)o), 

-sin(t-|-0o), 



(5) 
(6) 



(7) 

(8) 



where x, y are the Cartesian coordinates normalized to 
c/uji in the plane of the rotation of the vector E, t nor- 
malized to p is the electron momentum normal- 
ized to mcaQ — eEg/uJi, and (j)o is initial phase. We are 
interested only in particle dynamics in the momentum 
space so far as the electric field is homogeneous. We in- 
troduce Tp universally for all species of the particles as 
the angle between — E and p. Small and positive values 
of mean that particle lags behind the vector — E in 
the rotation. Then we have: px = —pcos{t + (jio — ip)i 
Py — —psin{t + ipQ — ip), and Eqs. ([7]) and ^ can be 
rewritten as follows 



dt 



1 - 



sin?/; 
P 



dp 



= cos ?/;. 



(9) 
(10) 



These equations can be reduced to the single equation in 
the plane p — ip: 



d(sin -ijj) 



= 1 - 



simp 



dp p 
The solution of the equation is the following: 



sin'0 = 



2 p \ 



sin xpQ 



Po 
2 



(11) 



(12) 



where po and "00 are the initial values of the electron mo- 
mentum and the angle, respectively. The electron trajec- 
tories ([T^ are shown in Fig. [T] (a). The positron trajec- 
tories (Fig. [1] (b)) can be obtained from Eq. ([T^ by the 



substitution ip — )■ ■0— tt. What this means is the positrons 
are accelerated in the opposite direction to the electron 
acceleration direction. The photon momentum does not 
change during the propagation, and the angle ip for a 
photon grows linearly in the course of time = ipo + t- 
Corresponding trajectories are shown in Fig. [l](c). 

The restoration of the particle energy in the cascade 
can be illustrated using Fig. [TJ Let us consider initially 
immobile electron. At first, from the point ip ^ 0, j = 
it is being accelerated along the bold line in Fig. [T] (a) 
(separatrix) : p ^ t and smip = p/2. Then it emits a 
photon and shifts left: its angle ip remains the same and 
its gamma-factor decreases due to the photon emission. 
Next according to the phase flow the electron comes again 
to the separatrix. This also explains why the electrons 
are concentrated near the separatrix. 

As we will see below, the electron distribution function 
drops along this line due to the following reasons. First, 
some of the electrons emit photons and leave the separa- 
trix. Second, the acceleration of the electrons that have 
higher energy starts earlier in the average than the accel- 
eration of the electrons with lower energies. The number 
of particles is increasing exponentially with time in the 
cascade, hence, the number of particles that gain higher 
energy is lower. 



III. BASIC EQUATIONS 

The cascade kinetic equations have been derived a long 
time ago to study development of showers initiated by 
cosmic rays in the Earth's atmosphere [H, [s^. The 
shower (avalanche) develops due to the bremsstrahlung 
and pair production at the particle interaction with nu- 
clei. The solution of the equations has been obtained in 
Ref. [13 . Later the solution method has been extended to 
include multiphoton processes [40| and to study the elec- 
tromagnetic showers in strong magnetic field [41t] . The 
kinetic approach have been used in numerical simulations 
to analyze the interaction of the intense laser pulse with 
the electron beam [l^. The kinetic equations have been 
also formulated for the problem of the cascade develop- 
ment in the laser field and have been used to derive ra- 
diation reaction force [13] in classical limit in Ref. [37| . 
However, up to now the kinetic cascade equations have 
been analytically solved only in one-dimensional cases. 
On the contrary, we consider theoretically and numer- 
ically the two-dimensional cascade. Moreover, we will 
take into account that the overall energy of the cascade 
particles is growing and can be much more than the seed 
particle energy. 

The cascade kinetic equations in arbitrary field can be 
written as follows: 
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Figure 1: The trajectories in p — t/j space for (a) electrons, (b) positrons and (c) photons moving in the rotating electric field. 
Photon emission and pair production are not taken into account. The bold lines correspond to the trajectories which separate 
bounded and open trajectories. 



9t/e(r,p) + ^V/e(r,p)+Vp{/e(r,p) xFe(r,p)}- / ^^/^(r, pXp' ^ p) dp' 
P ^p'IIp P 

+ / ^jj^fe{r,p'}w{p' ^p)dp'- Mr,p)w{p^p')dp\ (13) 
Jp'llp P Jp'Wp 



dtfp{r, p) + -V/p(r, p) + Vp{/p(r, p) x Fp(r, p)} 



,D-1 



— f^{r,p')w{p' p' -p)dp' 



+ / ^D^fpi^^ p'Mp' ^ P) dp' - / /p(r, p)u;(p -> p') dp', (14) 

Jp'llp P Jp'Wp 

9t/^(r,p) + -V/^(r,p) = - / /^(r,p)w(p-> p')dp' 
P Jp'Wp 

r- iD-l P ,D-1 

+ ^Tj^fe{r,p')w{p' ^p' -p)dp' + P^f^{r,p')w{p' ^p' -p)dp', (15) 

Jp'Wp P •'p'IIp P 



where r is the particle coordinate normalized to c/tui, F 
is the Lorentz force normalized to aomcuji, D is the space 
dimension, /e, fp and are the distribution functions of 
the electrons, the positrons and the photons, respectively, 
which are normalized such that 

y"/e,p,^(r,p)d^rd^p = iVe,p.7, (16) 
where Ng, Np, is the number of the electrons, the 



positrons and the photons, respectively. w{p' — >■ p)dp is 
the probability in time unit for the electron with momen- 
tum p' to emit a photon and to switch to the state with 
the value of momentum in the interval (p, p + dp) and the 
direction of the momentum parallel to p', w{p' — )■ p)dp 
is the probability in time unit for the photon with mo- 
mentum p' to decay with creation of the electron with 
momentum in the interval {p, p + dp) directed parallel to 
p' and the positron with momentum p' — p. Here we as- 
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sume that the photon emission and the pair production 
occur in synchrotron regime, and we neglect the angles 
less or about than I/7. Dependences of fe,p,-y, F on t and 
w, w onF are not written for the simplicity. 

The RHS of Eqs. (IT^ - p3|) describes photon emission 
and pair production. Namely, the first term in the RHS 
of Eq. (IT3I) characterizes creation of the electrons with 
momentum p by decay of the photons with momenta p' > 
p, the second and the third terms respectively describe 
the increase and decrease of the number of the electrons 
with momentum p due to photon emission. The first 
term in the RHS of Eq. (fT5|) characterizes photon decay, 
and the last two terms describe emission of new photons 
by the electrons and the positrons. 

The kinetic equations in the rotating electric field ([5]), 
([6]) in the plane geometry {D = 2) can be written as 
follows : 



dpx 



T sin(t + </.o)^%^ = / ^^(P')^(P' ^ P) dp'+ 

OPy Jp'llp P 

[ -fe,p{P>iP' '^P)dp' ~W{p)Up{p), (17) 
VIIp p 

+ [ P^/,(p')u;(p'^p'-p)dp' 
VlIP P 

+ [ ^MpW ^p' -p)dp\ (18) 
VlIP P 

where the relation ^(p' — > p' — p) = w(p' — > p) is taken 
into account, " in corresponds to the equation for 
the electron distribution function and " +" corresponds to 
the equation for the positron distribution function. 



/e,p,7(p) = J fe.pni'^i P) c^^r, 
W{p) = J w{p~^ p') dp', 

W{p) = / w(p ^ P')dp'. 



(19) 
(20) 
(21) 



In order to eliminate the dependence of the force and 
the probabilities on time in Eqs. p7)) and ([T5| . we make 
the change of variables Px,Py ^ p, ip as follows. Intro- 
ducing particle distribution functions in new variables 

ge,p,j{t,P,1p) =Pfe,p,j{t,Px,Py) = 

Pfe,p,i{t, -pcoa{t + 00 - V), -psin{t + 00 - V')), (22) 



we obtain: 



^9e-P-^ ^^fi^^P^l ^9e,P;^ 



dt 



dt 



(23) 



cosit + 00 ) 



— + sm{t + 00 Hr^ 
opx opy 

sin ipdge.p d ge,p , 

= — 2 — H-r-cos-07^ -. (24 

p-^ Oip op p 



Then, from Eqs. dTT]) and (US]) we have: 

dge.p _ dge^p _^ sinV^ gffe.p zppcOGi/-' ^ 9e,p{p,^) ^ 

dt dip p dip dp p 

9y{p', ip)w{p' p, ip)dp'+ 

{p', ^)wip' -> p, iP) dp' - Wge,p, (25) 



[ffe(p', V') + 9p{p\ ^)] w{p' ^p'-p, V) dp', (26) 



where ge,p,-y, W and W are defined at the point (p, ip). 

The probability rates w, w contained in Eqs. ([25]) . ([26]) 
can be found in the framework of the quasiclassical theory 
MM: 



w{p' 



a 1 



SldQ p' 



2r] 



2/3 



w{p' 



where 



(l_^)2/3 

a 1 



(l-^)^^^^x'^/^Ai'(.) 



7I/3 



, (27) 



£iao p'' 
X ( 2772/3(1 _^)2/3 



m)d^+ x''^'Ai'(i<) 
1 

" 771/3(1- j^) 1/3 



(28) 



1^ 

vx' 



'7(1 - v)x' 



2/3 



2/3 



(29) 



(30) 



V = p/p\ x' = SiOoP' sinip is the quantum parameter 
for the initial particle, e; = tuvi/mc^ is the frequency of 
the field rotation in Compton units, a — e'^ /he is the fine 
structure constant, Ai(^) and Ai (£) are the Airy function 
and its derivative, respectively [43]. Eqs. (|25)) -p8 l) are 
the complete set of equations which describes evolution 
of the distribution functions in the momentum space for 
the electromagnetic cascade in the rotating electric field. 



IV. SCALABILITY OF THE CASCADE 
EQUATIONS 

Eqs. and (pS)) governing cascade dynamics depend 
on two parameters: period of the electric field rotation 
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Figure 2: The distribution functions of (a) the electrons and 
(b) the photons of the electromagnetic cascade developing 
in the rotating electric field with the normalized amplitude 
1.6 X 10^ and A; = 0.8 fim. 



ao 



and the field strength. However, it can be shown that 
if Eq 3> aEcr the number of the parameters can be 
reduced. It is shown in Ref. [s^l that if Eq ^ aE^ 
and A; > 2'Kh/mca^ ~ 50 nm the characteristic time be- 
tween consecutive photon emissions by the cascade elec- 
tron (positron) are much smaller than the period of the 
electric field rotation 27r/aj;. In other words, the number 
of the cascade particles increases in many times during 
the period of the electric field rotation, F ^ 1, where F 
is the cascade growth rate normalized to w;. The number 
of cascade particles increases exponentially N (x exp(Ft). 

The simplification of the kinetic equations and the re- 
duction of the number of the governing parameters are 
possible due to the fact that in this case the cascade par- 
ticles are concentrated in the regions Isin?/!! <C 1. As 
we will see later it follows from numerical simulations 
that the condition | sin -01 <C 1 holds for the most of the 
particles produced in the cascade. This can be qualita- 
tively explained as follows. Eq. (fTO| implies that p ^ 1 
achieves when the electrons or positrons are accelerated 
within the time interval about the period of the electric 
field rotation. If Eq ^ ctEcr, the electron emits a lot of 
photons within this time interval and loses a great part of 
the gained energy. Hence, the electron does not gain high 



energy from the field and its characteristic momentum is 
small p <C 1. In accordance with Eq. ([9]), the angle ip 
between the electron velocity and the direction opposite 
to the electric field direction, changes very quickly until 
I sin S P ^ 1 ■ The photons emitted by the electrons 
mostly have | sin -01 <^ 1. Furthermore, for the photons 
dip/dt = 1, so they decay before the angle ip changes 
significantly. Therefore, it is natural to assume that the 
distribution functions of the electrons, the positrons and 
the photons are substantially nonzero only in the region 
I sin -01 < 1. 

Using the fact that Isin?/)! <^ 1, the kinetic equations 
can be simplified as follows. |sin0| ^ 1 if « or 
-0 « TT. Let us, for the specificity, consider the region 
^ 1. In this region the following substitutions can be 
made in Eq. (j25|) : cosip — 1, simp — -0. In the resulting 
equation as well as in Eq. (j26p (with probabilities from 
Eqs. ^7^1, ^S^) the following substitution 



p = Pp, 



(31) 
(32) 

(33) 



does not change anything except parameter x ~ 
Sio^ptp — eiUQ^'^ptp, where j3 is arbitrary quantity. 
Choosing 



1 



1/2 ■ 



(34) 



yields x = ^"0- The resulting equations for the distribu- 
tion functions depend only on the parameter 



1 Eq 

a Ecr 



(35) 



They do not depend on the frequency uji and, hence, 
on the corresponding wavelength A; = 2ttc/uji. Thus, 
the dependence of various cascade characteristics on A; 
is determined only by the scaling Eqs. (PT|) - (IM1) . For 
example, it follows from Eq. ((3T|) that for the certain 
field amplitude and, hence, the intensity / — cEq/Att, 
the growth rates at two different wavelengths relates as 
follows: 



r(j,A 



r(/,A,,2) 




(36) 



This formula is in good agreement with the results of the 
numerical simulations as shown in Sec. VII. 



V. SCALABILITY OF THE KINETIC 
EQUATIONS FOR PARTICLES WITH x > 1 

Further simplification can be made if we consider high- 
energy "tails" of the distribution functions with x ^ 1- 
This let us to simplify the formulae for the spectral dis- 
tributions of the probabilities Eq. (gT]), Fist, in the 
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2lg{///25) 



Figure 3: The growth rate of the electromagnetic cascade in 
the rotating electric field as a function of J = cEq/Atv and 
A; = 2nc/uui. Here Eo is the amplitude of the electric field, uji 
if the cyclic frequency of the field rotation, I25 = fO^^ W/cm^. 



region where the probabilities are substantially nonzero, 
the derivative and the integral of the Airy function are 
of the same order: /jAi(^)d^ - Ai'(x) if < < 1. 
Second, 



27]' 



,2/3 



771/3(1 _^)l/3 



(l-7?)2/3 ' ^1/3 J 



-2r;2/3(l_r;)2/3 



2, (37) 



0.8, (38) 



if ?7 G (0,1). Hence, as soon as x ^ 1; the terms in 
Eqs. (p7| and (p8)) that are proportional to the integrals 
of the Airy function can be neglected. Furthermore, >c is 
much smaller than unity on a wide interval 1 /x' ^ ^7 < 1 , 
and k 1 ii 1/x' ^ rj and 1 — 77 ^ Then in 

Eqs. ((27|) and ((28)) we can set k = >c = 0. However, this 
lead to appearing of the singularities in the expressions 
for w and w, namely w oc l/ij^^^ at ry 0, oc l/rj^^^ 
if r/ ^ and w cx 1/(1 — ij)^/"^ if r/ 1. Nevertheless, 
these singularities are integrable and give contributions 
to the total probability rates that are about x'^^^ ^ 1 
times smaller than the total probability rates. 
Finally, for x' ^ 1 and ^ 1 we have: 



w(P P) 



Ai'(O) 



aa. 



1 



w{-p' p) = _ Ai'(O) 



gl/3 p/4/3 ^1/3(1 _^)2/3' 

aal'^ 7/.2/3 + (1 - 7]f 

p/4/3 ^1/3(1 _^)l/3' 

e/ p' 



(39) 
(40) 
(41) 



and the same expression for W except substitution v for 
V, where 



Ai'(O) 



■dr]K, 1.46, 



(42) 



/o ,71/3(1 -,y)2/3 

Ai^(0)/ \v3^i%v'3 ^^^0-38. (43) 



Furthermore, Eq. ()26p can be used in the same form and 
Eq. ([^5)) can be simplified as follows: 

dge,p dge,p . V' dge,p d ge,p , 

— — — ± — 

dt dip p dip dp p 

gj{p\ ip)w{p' ^ p, -0) dp'+ 

{p',^)wip' ^p,ij)dp' -Wge^p. (44) 

Similar equations can be written in the region x 3> 1, 
I?/)- 7r| < 1. 

The substitution: 



-1/4 



i3/4ay 

.1/4 



1/4^' 



a3/4ai/4 



(45) 
(46) 
(47) 



in Eqs. (^5]). ([23)- (El) yields the following parameter- free 
equations: 

dge,p dge,p , V' dge,p _ 9 5e,p , 

=^ = =^ ± — =^ ^ p — ^ — + 

dt dip p d'4) dp p 

g^{p', Tp)w^f{p' p, Tp) dp'+ 
ip', i>)wip' ^ p, ^) dp' - Wge,p, (48) 



P 
00 



dg^y dg^ - 

[ge{p', ^) + V^)] -^P' -P. ^) dp', (49) 



ip 
where 



w(p' ^ p) = - Ai'(O) 
w^(p' ^ p) = 



j;2/3 1 + 



p/4/3 ^1/3(1 _^)2/3' 
2 



^ ^p'4/3 ,,1/3(1 _^)l/3' 



(50) 
(51) 
(52) 



W^^(p,V)=^)|j7j, (53) 
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Figure 4: The dependence of the growth rate on A;: the nu- 
merical data (dots) and the fits T/uji = Ay/Xi for the data 
(lines). The amplitude A is chosen such that the RMS de- 
viation from the numerical data is minimal. The error bars 
show the dispersion of the numerical data for 7 runs of the 
simulation. Line 1 corresponds to the data obtained for 
/ = lO^'* W/cm^, line 2 corresponds to 7 = 10^^ W/cm^ and 
line 3 corresponds to J = 10^^ W/cm^. 




Figure 5: The dependence of the growth rate on /i = 
Eo/aEcr- The numerical results are presented by the solid 
black dots. The error bars show the dispersion of the numeri- 
cal data for 7 runs of the simulation. The solid lines represent 
fits r = Afi^^'^. The amplitude A is chosen such that the RMS 
deviation from the numerical data is minimal. Line 1 corre- 
sponds to numerical data for A; = 10 nm (squares), line 2 
corresponds to A; = 5 /xm (triangles) and line 3 corresponds 
to A; = 1 nm (circles). 



The obtained equations can be used for further analysis of 
the particle distribution during the cascade development. 
However, these equations do not completely determine 
the cascade dynamics, because they describe only the 
particles with high x- For example, the particles with 
X 1 can also breed and give some contribution to the 
growth rate. 

Let us assume that the growth rate is determined only 
by the particles with x ^ 1, IV'I 1 f^nd It/" — 7r| <C 1. In 
this case the dependence of F on the electric field strength 
and the wavelength can be determined from Eq. (j45l) as 
follows: 



r(/xi, A;,i) 
r(M2, A;,2) 



,1/4x1/2' 



(54) 



VI. ANALYTICAL SOLUTIONS AND ENERGY 
SPECTRA 



Here we again restrict ourself to the limit Eq S> aEcr, 
and consider the "tails" of the distribution functions that 
correspond to x ^ 1 and, additionally, to (p) in or- 
der to estimate the terms in Eqs. (|^5)) . and obtain 
approximate solutions for the stationary particle distri- 
bution functions as well as for the energy spectra of the 
electrons, the positrons and the photons. Here (p) is 
the average particle momentum and the stationary dis- 
tribution means that the distribution function shape is 
approximately conserved during the evolution, but the 
number of particles grows exponentially: 



5e,p,7 oc exp(rt). 



dt 



7- 



(55) 
(56) 



where /i = Eo/aEcr- The dependence Eq. ([S3]) and 
the scaling Eqs. p5)) -(|47 | was also derived in Ref. ^] 
via simple estimations and dimensional analysis of the 
kinetic equations. However, the numerical simulations 
show that in a wide range of the parameters Eq. (I54 p is 
not valid (see Sec. VH and right plot in Fig. 7 in Ref. |37| ) 
and can be used only for estimations with accuracy about 
50%. Thus, the growth rate F is not completely deter- 
mined by Eqs. p5|) - ([55)) that describe only particles with 
X ^ 1. However, the growth rate contributes to these 
equations through time derivatives. Hence, for further 
analysis of Eqs. (P5|) - P5|) we assume that F is known (for 
example, from the numerical simulations). 



The numerical simulations shows that the distribution 
functions approach such dependence in time interval 
about some 1/F after the cascade starts. 

The main source of the electrons with {tpl ^ 1, p ^ (p) 
is their transfer from the region p < (p) due to the accel- 
eration. So, electrons are concentrated along the separa- 
trix p/2 owing to the phase flow. We will see that the 
'0-width of the electron distribution function at p 3> (p) 
is much smaller than the angular width of the photon 
distribution function and much smaller than p. Because 
of that we derive the electron energy spectrum and do 
not derive the angular shape of the electron distribution 
function. 
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Furthermore, we will show that particle spectra 



9e 



(57) 



decrease exponentially with the increase of p due to trans- 
fer of the electrons with high momenta to the region of 
lower momenta because of frequent emission of the pho- 
tons. Besides, the characteristic scale of the decrease is 
less or about the average momentum (p). Hence, the 
electrons that acquire momentum p due to the photon 
emission mostly have momentum p' before the emission 
in the following narrow interval: p < p' ^ P + (p)- The 
electrons that initially have momentum p are distributed 
on wide interval p' £ (0,p) with the average momentum 
about p after the photon emission, because the character- 
istic p'-scale of the function w{p — )■ p') is the same order 
as p. The last two terms in Eq. (^5)) describe the photon 
emission. The last but one term in Eq. (|48l) describes the 
inflow of the electrons from the region of high momenta 
p' > p, and the last term describes the outflow of the 
electrons to the region of momenta less than p. Sum- 
marizing, we obtain that the last but one term is about 
p/{p) 3> 1 times smaller than the last term and can be 
described in the framework of the perturbation theory. 
It can be shown by similar reasoning that the term with 
in Eq. can be neglected. 

Integrating Eq. for the electrons over the interval 
^ S {p/2 — A'!/;,p/2 -|- A-i/i), where Aip <^ p/2 and much 
higher than the T/^-width of the electron distribution func- 
tion, and neglecting g^, we obtain: 



dhe 
dp 



-TK 



vp 



1/3 



22/3 

2i/3Ai'(0)pi/3 



he{p[l+^]) 
^2/3 



(58) 



where we simplify the formula for w{p' — ^ p, V") taking 
into account that ip p' /2 and p' — p ^ (p) <^ p. As 
it was already mentioned, the characteristic scale of the 
integrand in the last term of Eq. (|58|) is ^ ~ (p) /p <^ 1 
and this term can be treated as perturbation. In the 
zeroth order of the perturbation theory we have: 



he{p) = ho exp 



-Tp 



3vp' 



,4/3 



28/3 



(59) 



where ho is a constant. Then, using the following ap- 
proximation derived by the expansion of the terms in the 
exponent in Tailor series 



/le(p[l +C]) ~ /le(p)exp 



vp 



;4/3 



22/3 



we obtain: 

1 dhe 
he dp 



vp 



1/3 2l/3Ai'(0)6_2/3 



22/3 



[f -f wpl/3/22/3] 



1/3 ' 



(60) 



(61) 



where 



6,; - 



/ rexp(-OrfC, 
Jo 



(62) 



and 5_2/3 ~ 2.68. Finally, from Eq. (|6ip we have: 



, j ^_ 3«p-4/3 9Ai'(0)b,2/3f«/\ , 
he(p) = ftoexp < -Ip TTTTT^ = „ X 



28/3 



5 X 22/3^3 



2/3 



(5p2/3-6pi/3 + 9) -9 



, (63) 



where p — v^p/AT^. 

The positrons are decelerated in the region ^ 1, 
so, the only source of the photons in the region \^p\ <C 1, 
p ^ (p) are the electrons, hence, we can neglect gp and 
assume p' — p ^ p' in Eq. P^j) . that yield: 



vj;^/^ _ 2^/3 Ai'{0)he{2{l;) 



pi/3 {2x1; - p) 



1/3 



(64) 



In order to obtain this equation we also use the electron 
distribution function in the following form: 



ge{p'.i^)^he{p')5{{p-p'l2), 



(65) 



where 5{x) is the Dirac delta function. Obviously, the 
photons, emitted by the electrons in the region <C 1, 
are situated in the p — ip space above the line xp = p/2. 
As far as we assume that the photons decay before they 
pass the distance about tt along i/j-direction, we can set 
gjiPiP/"^) = 0. Additionally, making of use the following 
approximations in the solution of Eq. (|64p 



he{2tp) K, he{p) exp 

*--(! 

we obtain for p/2 < ■0 <C 1 

9j{p,i') = 



r 

5/3 



vp 



1/3 



22/3 
Jp2/3 

3 X 



(2^-p) 

22/3- ^ f 



(66) 
(67) 



Ai'(0)/le(p) 



P 



1/3 [f - 



(2l/3t;-w/22/3)pl/3 

^ exp(-0 



X exp 



(*-f) 



^1/3 

3?) 



^5/3 



p5/3 
W3 



(68) 



where 



+ (21/^1; - 5/2^/3) pi/3] (^ _ p/2) . (69) 



Further simplifications can be made if we assume that 
2i/3^pi/3 ^ f and take into account that v/v 0.26 < 
1. In this case the characteristic scale of the integral in 
Eq. is much smaller than the scale of the exponent in 
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Figure 6: The electron (line 1) and the photon (line 2) dis- 
tribution functions obtained by the numerical simulation of 
the electromagnetic cascade developing in the rotating elec- 
tromagnetic field with ao = 8 x 10^ and A; = 0.8 fim for 7: 
7/ao — 0.05 (a) and 7/ao — 0.15 (b). Lines 3 correspond 
to the analytical expression for the photon distribution func- 
tion Eq. (|68|) . All plotted functions are normalized to their 
maximal values. 



Eq. ([55)1 . Thus, integrating Eq. over ip, we assume 
that the integral in this equation is equal to its maximal 
value 6-1/3 if > 0. We also suppose that the positron 
distribution function concentrated generally near ^ = tt 
gives the same photon distribution as Eq. ([55)) except 
the substitution -0 — tt for ip. Furthermore, in order to 
advance simplification we can use Taylor's expansion in 
the exponent of Eq. (1551) similarly to Eqs. (p5t and (|57|. 
Finally we obtain the following expression for the photon 
energy spectrum: 




Figure 7: The electron (line 1) and the photon (line 2) dis- 
tribution functions at 7/ao =0.1 obtained by the numeri- 
cal simulations of the electromagnetic cascade developing in 
the rotating electromagnetic field with A; = 0.8 fim and (a) 
ao = 3.2 X 10'', (b) ao = 1.6 x 10^ (c) ao = 8 x 10^ Lines 3 
correspond to the analytical expression for the photon distri- 
bution function Eq. (|68p . All plotted functions are normalized 
to their maximal values. 



2Ai'(0)6_i/3/le(p) 



pl/3 [f -f (2l/3„-_i)/22/3)pl/3 

1 



2/3 



[f -^{ipl/3/22/31 



(70) 



VII. NUMERICAL SIMULATIONS 

In this section we compare the analytical expressions 
for the electron energy spectrum ([63)) . the photon en- 
ergy spectrum ()70p and for the photon distribution func- 
tion (I68|) with the results obtained using the numerical 
model described in Ref. [28]. Furthermore, using the nu- 



merical values of r, we compare various terms in these 
expressions and simplify them for a wide range of param- 
eters. 

The electron and the photon distribution functions in 
electromagnetic cascade are shown in Fig. |2j These func- 
tions are obtained in the numerical simulation performed 
for the field parameters oq = 1.6 x 10^, A/ = 0.8 /im. The 
numerical simulation gives T/uji w 12.1 for these parame- 
ters. The shape of the distribution functions is typical for 
the cascades with T 1 . The electron distribution func- 
tion at p 3> (p) is concentrated near the line ip = p/2, 
and the photon distribution function is nonzero upper 
this line, that agrees well with the qualitative descrip- 
tion given in Sec. VI. Furthermore, it follows from Fig. |5| 
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Figure 8: The electron (upper lines) and the photon (lower 
lines) energy spectra for the electromagnetic cascade devel- 
oping in the rotating electric field with A; — 0.8 fim and (a) 
ao = 3.2 X 10'', (b) ao = 1.6 x 10^ (c) ao = 8 x 10^ Solid 
lines correspond to the numerical results and the dash-dotted 
lines correspond to the analytical expressions H63p and (|70|) . 



and the results of the numerical simulation that almost 
all particles are located in the regions p <^ 1, ^ 1 
and p < 1, IV- - 7r| < 1. 

The dependence of the growth rate F on the electric 
field parameters is shown on Fig. [3l In every node of the 
surface F is calculated 7 times and the result is averaged 
out. It is seen that F grows with the increase of the 
intensity and the wavelength, and the growth of F due to 
the wavelength increase is much steeper than the growth 
due to the intensity increase. Therefore, it is preferable 
to use longer wavelengths in the experiments devoting to 
the cascades observation. Besides, the limitation of the 
attainable intensity due to prolific pair production is less 
strict for shorter wavelengths. 

In Fig. U the dependence of the growth rate on the 
wavelength for three different values of the intensity is 



shown. The rows of dots show the results of the numeri- 
cal simulations. The error bars show the dispersion of F 
obtained for 7 simulations with the same initial parame- 
ters. The cascades in the simulations are initiated by a 
single immobile e"'"e~ pair. Thus, at initial stage of the 
cascade development when the number of cascade par- 
ticles is small the shape of the distribution functions is 
being modified. Because of this the growth rate fluctu- 
ates significantly at the initial stage. We measure the 
growth rate as F = d{\nNe)/dt at the later stage of 
the cascade development when the number of particles 
is high and the growth rate becomes constant. The dis- 
persion of F, obtained in the simulations, is explained by 
the finiteness of the particle number used in the simu- 
lations. The solid lines in Fig. 2] are the fits described 
by the equation F = A\/Jq, where A is chosen for ev- 
ery row of the numerical data such that the RMS dif- 
ference between the numerical data and the fit is mini- 
mal. The line 1 corresponds to the numerical data ob- 
tained for / = 10^^ W/cm^, the line 2 corresponds to 
the numerical data obtained for / = 10^^ W/cm^ and 
the line 3 corresponds to the numerical data obtained for 
/ = 10^^ W/cm^. It is seen from Fig. |3]that the theo- 
retically obtained dependence of the growth rate on the 
wavelength Eq. (PSI) is in good agreement with the results 
of the numerical simulations. 

The dependence of the growth rate on the electric field 
strength is shown in Fig. [S] Black points present the 
results of the numerical simulations, and the solid lines 
are the RMS fits F = Afi^/"^ that correspond to Eq. ((54)) . 
The line 1 corresponds to the squares that illustrate the 
numerical results obtained for A/ = 10 /im, the line 2 
belongs to the triangles that illustrate the numerical data 
for Xi — 5 /im and the line 3 corresponds to the numerical 
data obtained for Xi — 1 fj.m (circles). It follows from 
Fig. [5] that Eq. ([54]) can be used only for estimations 
of F with accuracy about 50%. However, Eq. is 
derived under the assumption that x ^ 1 for the most 
of cascade particles, but actually many cascade particles 
can have x ^ 1- These particles can give substantial 
contribution to the cascade dynamics and, hence, to the 
cascade growth rate. Indeed, numerical simulations show 
that for / = 10^^ W/cm^, A; = 1 /im the portion of the 
electrons with x < 5 is about 30% and this portion for 
/ ^ 10^5 W/cm2, A; = 1 /im is about 90%. 

The comparison between Eq. (|68p and 5ph(V') at cer- 
tain value of p obtained by the numerical simulations is 
given in Fig. |6l The line 1 and the line 2 correspond to 
the electron and the photon distribution functions, re- 
spectively. The line 3 corresponds to the solution (|68l) . 
The distribution functions are normalized to their max- 
imal values. The parameters of the electric field are: 
ao — 8 X 10^, A/ = 0.8 /xm. The distribution functions 
gc-yiipjp) are shown at p = 0.05 (plate a) and at p = 0.15 
(plate b) . The growth rate used in Eq. ([68)) is calculated 
numerically: F/o;; — 24 A. It is seen that the electron 
distribution function is much narrower than that of the 
photons. Furthermore, the numerical results are in good 
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agreement with Eq. (j68l) . Similar cuts of the distribution 
functions are shown in Fig. [T] The panel (a) corresponds 
to flQ — 3.2 X 10* and F = 6.29 obtained in the simulation, 
the panel (b) corresponds to oo = 1.6 x 10^ and F — 12.1, 
the panel (c) corresponds to oq = 8 x 10^ and F = 24.1. It 
is seen from the plates (b) and (c) that Eq. (|68l) is in good 
agreement with the results of the numerical simulation. 
However, the photon distribution function obtained the- 
oretically reaches its maximum earlier than the photon 
distribution function obtained in numerical simulations. 
Possible explanation of this small discrepancy is the fol- 
lowing. In order to obtain Eq. (1551) we assume that the 
electron distribution function depends on the angle '0 as 
follows; ge{p,4') cx 6{ip — p/2), however, the numerical 
simulations show that ijj-width of the electron distribu- 
tion function is finite and big amount of electrons have 
"0 slightly higher than p/2, hence, the maximum of the 
photon distribution function is also slightly shifted to the 
region of higher angles. Besides, for the plate (a) the con- 
dition p ^ (p) is not fulfilled, the shape of electron distri- 
bution function in this case becomes more complex than 
in the case p (p) and Eq. (|68|) gives rather qualitative 
approximation for the photon distribution function. 

In Fig. [8] the electron and the photon energy spectra 
are shown. The solid lines correspond to the results of 
the numerical simulations and the dash-dotted lines cor- 
respond to Eqs. ([55]) and ([70]) . The top lines correspond 
to the electron distribution functions and the bottom 
lines correspond to the photon distribution functions. In 
Eqs. ((63|) and ([70|) we use F from numerical simulations. 
A; = 0.8 fim, plate (a): ao = 3.2 x 10*, F = 6.29, plate 
(b): ao = 1.6 x 10^ F = 12.1, plate (c): oq = 8 x 10^ 
F = 24.1. It is seen that the obtained analytical results 
are in good agreement with the results of the numerical 
simulations. 




Figure 9: The ratio of the growth rate, F, to the fit Eq. (ffT|) 
(dots) and unity level (line). F is computed for the intensity 
range / — 10^^ — 10^^ W/cm^ and the range of wavelengths 
A = 1 — 10 /im with the step 1 fim. 



We can propose the following fit for the growth rate: 



T/uji ^ 0.65 X Xl^^lfim 



Jo / Jo 



(71) 



where /q = 6 x 10^^ W/cm^, the magnitude and Iq was 
chosen in order to obtain good agreement with F from 
the numerical simulations. The ratio of F obtained in 
the numerical simulations of Fig. [3] to the fit Eq. ([7T|) 
(circles) is shown in Fig. [HI It is seen that Eq. (|71l) has a 
satisfactory accuracy. 

Making of use Eq. ((7T|) we can estimate the terms in 
Eqs. dSS]), dMl) and dZO]). First, for f we obtain the fob 
lowing approximation: 



VIII. FURTHER SIMPLIFICATIONS 

The obtained analytical expressions for the particle 
distribution functions are quite complex. Besides, these 
expressions contain the growth rate F that is not derived 
in the framework of the presented theory. In this section 
we compare various terms in the obtained expressions 
and derive asymptotic (p — >■ c») expressions for the par- 
ticle spectra that are more simple than Eqs. (|63| . ([70| . 
In order to obtain these asymptotic spectra we keep in 
Eqs. and ([7(7| only the terms that describe the slope 
of the logarithmic spectra and the ratio between the elec- 
tron and the photon energy spectra. The slope and the 
ratio can be probably measured in the future cascade ex- 
periments that allows estimating of the field magnitude 
in the cascade region and the growth rate. Furthermore, 
we find the fit for F guided by the results of the nu- 
merical simulations in a wide range of the parameters 
(A, = 1 - 10 fim, I = 10^5 - 10^9 W/cm2). 



0.22 Ig 



lo' 



(72) 



Thus, the ratio of the first and the second terms in the 
exponent in Eq. is 



28/3fp ^ lg///o ■ 



(73) 



For the considered region of parameters \gl / Iq « 1 — 5. 
Therefore, asymptotically at high values of p this ratio is 
much more than unity and we can neglect the first term 
in the exponent in Eq. (|63p . Besides, the last term in 
the exponent in Eq. ([63l) relates to the second term as 
follows. 



0.3 



F4/3p4/3 



(pl/=^ + l) 



2/3 



X (5p2/3 - 6pi/^ + 9) - 9] , (74) 



where p is introduced right after Eq. (j63|) . This expres- 
sion decreases with the increase of p « Q.%p/T^ , hence. 
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neglect 2 2/3^^1/3 ^ 0.24p^/'^ in comparison with T. The 
resuhing asymptotic form of the photon spectrum is 



Figure 10: The electron (upper solid lines) and the photon 
(lower solid lines) energy spectra obtained by the numerical 
simulations of the electromagnetic cascade developed in the 
field characterized by A = 0.8 /xm and ao = 3.2 x lO"* (plate 
(a)), ao = 1.6 X 10^ (plate (b)), ao = 8 x lO'' (plate c). The 
corresponding asymptotic expressions (|75p and (|77|l are shown 
by the dash-dotted lines. 



for sufficiently large values of p the last term can be ne- 
glected for simplification. Therefore, the electron spec- 
trum Eq. can be reduced to the following asymptotic 
form: 



— ha exp 



(75) 



where p = a ' Oq 7 I.37/7* and 

7, = 2.4 x 10~^Ay^[/ini]/^/*^[W/cm2]. (76) 

In order to simplify Eq. ([70)) we neglect 5/2 sa 0.19 
in comparison with v sa 1.46 and we neglect T in com- 
parison with 2^/'^vp^^^ w 1.8p^^^. However, we does not 



0.47/le,s(7) 



p5/9 [f + 0.24pi/3] ' 
where 

-j^ 1/4 -3/4 -l/4-p / 

323 



(77) 



A;^/^[^m]Ji/8[w/cm2] 



(78) 



The comparison of Eqs. (|75l) and (1771) with the results of 
the numerical simulations is given in Fig. 1101 



IX. SUMMARY AND DISCUSSION 

In this paper we study the development of the electro- 
magnetic cascade in the rotating electromagnetic field. 
First we derive the kinetic equations for the electrons, 
the positrons and the photons in the cascade. Second, 
we present the dimensional analysis and scaling for the 
kinetic cascade equations. In the limit when the growth 
rate is quite large F 3> 1 the angle between the momenta 
of the cascade particles and the electric field vector is 
small. This allows us to find how the cascade parameters 
such as the average particle momentum and the cascade 
growth rate depend on the frequency of the field rota- 
tion. It is shown that lower frequencies is preferable in 
the experiments devoting to the cascades observation. At 
the same time, the limitation of the attainable intensity 
due to prolific pair production is less strict for higher fre- 
quencies. Then, under the assumption that the most of 
particles have x ^ 1 the dependence of the cascade pa- 
rameters on the electric field strength is found. Besides, 
for the distribution function "tails" corresponding to the 
particles with high values of quantum parameter x ^ 1 
and with Lorentz factor 7 3> (7), the analytical solution 
is found. The obtained formulae for the distribution func- 
tions and for the energy spectra agree fairly good with 
the results of the numerical simulations. Furthermore, 
the fit for the growth rate based on the results of the 
numerical simulations in a wide range of the parameters 
is proposed. 

The electromagnetic cascades can arise in various real- 
istic field configurations, particularly, in the field of cir- 
cularly or linearly polarized standing wave (28j . In the 
field of the circularly polarized standing wave the cascade 
develops mostly in the _B-node, where i? ss 0. Near this 
plane the field configuration is very close to the homo- 
geneous rotating electric field. Certainly, some particles 
can leave the plane B = 0. However, we assume that 
this fact does not affect cascade dynamics much and the 
proposed theory fits well for the description of the elec- 
tromagnetic cascades in the circularly polarized standing 
wave. The estimations also demonstrate that the cascade 
can arise in the crossed time-dependent electric and mag- 
netic fields if the magnitude of the electric field is greater 
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than the magnitude of the magnetic field. In this case 
particle momenta are also concentrated along the distin- 
guished direction and some simplifications of the kinetic 
equations are also possible. Thus, the proposed model 
for the electromagnetic cascades in rotating electric field 
can be used for analysis of the cascade dynamics in some 
other field configurations. 

It follows from the numerical simulations that in the 
considered region of the parameters the big portion of 
the particles has x ~ 1 and can strongly affect the cas- 
cade development and the cascade growth rate. As far as 
the obtained distribution functions for the particles with 
X ^ 1 depend on the cascade growth rate, the proposed 
description of the "tails" of the distribution functions is 
not completely self-consistent. The contribution of the 
particles with x ^ 1 in F is especially substantial for the 
near-threshold cascades. Thus, the investigation of the 
distribution function component with x ~ 1 is important 
for experimental observation electromagnetic cascades in 
laser field at the lowest possible intensity level. 

However, it is important to note that the predictions 



of the proposed theory for the dependence of the cascade 
parameters on the frequency of the field rotation and for 
the distribution functions and the energy spectra of the 
high-energy cascade particles are in fairly good agree- 
ment with the results of the numerical simulations. The 
proposed model may be useful for the development of 
plasma sources and for the better understanding of 
some astrophysical processes. Furthermore, the obtained 
expressions might be also used to estimate the cascade 
parameters in the experiments by measuring the particle 
energy spectra. 
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